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Vlt«  optimality  (although  no  counter  example*  ara  known).  The  problem  la  ap- 
proached aa  one  of  decompoalng  a network  aubjact  to  conatralnta.  The  problem 
may  be  atated  aa: 
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atralnta, a mean-tlme-to-repalr  conatralnt  and  aatabllah  aparea  to  meat  an 
availability  conatralnt. 


A computer  routine  haa  been  developed  and  explained. 
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EVALUATION 

This  study  was  one  of  many  responsive  to  Technical  Planning  Objective 
«13  (Reliability)  performed  by  Syracuse  University  under  Contract  No. 
F30602-71-C-0312.  The  objective  of  this  study  was  the  development  of 
mathematical  models  for  selecting  a set  of  modules  which  represents  the 
division  of  an  electronic  system  resulting  in  the  lowest  life  cycle  costs. 
The  study  provided  a method  based  on  network  theory  which,  while  not  proven 
to  provide  the  optimal  solution,  will  certainly  provide  an  efficient  means 
for  reducing  system  life  cycle  costs  through  better  partitioning  in  design. 

Future  work  related  to  this  effort  will  be  the  development  of  other 
tools  for  lowering  life  cycle  costs  through  more  cost-ef fective  design,  and 
the  incorporation  of  these  into  Air  Force  standards  and  handbooks  to  assure 
their  application. 


ANTHONY  COPPOLA 
Project  Engineer 
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I.  INTRODUCTION 

The  objective  of  this  research  was  to  develop  a method  of  estab- 
lishing modules  for  equipments  in  a way  that  would  "minimize"  the  life 
cycle  cost  of  the  system.  The  construction  of  the  modules  was  to  be  as 
non-restrictive  as  possible  without  an  exhaustive  enumeration  of  all 
possible  modules.  (Total  or  exhaustive  enumeration,  while  feasible  for 
small  networks  would  be  infeasible  for  large  networks.  A network  of  ten 
nodes  may  contain  more  than  115,000  different  modules.)  For  all  practical 
purposes,  the  number  of  possible  modules  in  a 100  node  network  is  infinite 
Jensen  (1,2),  Caponecchi  and  Jensen  (3)  and  Caponecchi  (4)  have  proposed 
the  use  of  proper,  restricted  cuts  as  a basis  for  establishing  modules. 

In  this  work,  those  restrictions  have  been  removed.  Proper,  restricted 
cuts  are  too  restrictive  and  could  result  in  the  construction  of  modules 
that  are  extremely  non-optimal. 

The  solution  procedure  developed  in  this  research  requires  that  the 
network  be  partitioned  N times,  where  N is  the  number  of  nodes  in  the  net- 
work. The  "best"  of  the  N modular  arrangements  is  chosen  as  the  desired 
modular  structure.  The  "best"  structure  is  that  structure  which  has  the 
lowest  life  cycle  cost.  A computer  routine  and  numerical  example  are 
given. 

II.  NETWORK  STRUCTURE 

The  systems  under  consideration  must  be  capable  of  being  structured 
as  networks.  (Electronic  networks  are  prime  examples.)  Each  node  in  the 
network  must  consist  of  an  element  that  cannot  feasibly  be  decomposed. 

In  other  words,  each  node  can  be  called  a "smallest  functional  element"; 
that  is  an  SFE.  The  arcs  in  the  network  represent  the  interconnections 
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between  nodes  or  SFE's.  In  essence,  they  (the  arcs)  are  distinct, 
discreet  channels  by  which  the  nodes  communicate.  An  example  network 
is  shown  in  Figure  1. 

Networks  can  be  classified  in  several  ways.  They  are  graphs.  In 
fact,  they  are  geometric  graphs  and  may  be  planar  or  non-planar,  dir- 
ected or  non-directed  and  feed-forward  or  feedback.  The  developments 
that  follow  apply  to  any  network  that  fits  into  the  above  description. 
(For  a more  detailed  description  of  graphs,  see  Busacker  and  Saaty, 
Finite  Graphs  and  Networks,  McGraw-Hill,  1965,  Chapters  1 and  2.) 

III.  NETWORK  PARTITIONING 


Networks  can  be  divided  into  subnetworks  or  modules  by  two  basic 
types  of  partitioning.  The  two  types  are  1)  those  minimum  partitions 
that  disrupt  all  communcation  paths  from  source- to- sink,  and  2)  those 
partitions  that  do  not  disrupt  all  communication  paths  from  source-to- 
sink.  The  network  in  Figure  2a  is  partitioned  to  discrupt  all  paths 
while  the  same  network  shown  in  Figure  2b  is  partitioned  so  that  not  all 
paths  are  disrupted. 

The  first  type  of  partitioning  is  called  partitioning  by  proper 
cuts.  When  the  network  has  directed  arcs,  it  is  possible  to  use 

•» 

restricted  proper  cuts.  The  method  being  proposed  permits  any  combination 
of  botfh  types  of  partitioning.  The  partitioning  may  result  in  disjoint 
elements  (those  with  no  direct  interconnections)  being  assigned  to  a 


module. 


An  example  of  such  partitioning  is  shown  in  Figure  3. 
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labels  are  number  of  Interconnections  between  nodes 


Sff** 


IV.  PROPER  CUT  PARTITIONING 
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Given  a linear  graph  (network)  a proper  cut  is  defined  as  a set 
of  arcs  which  when  removed  from  the  graph  divide  it  into  exactly  two 

N 

connected  subgraphs.  For  a graph  of  N nodes  there  are  a maximum  of  2 

different  partitions  of  the  nodes  into  two  subsets  X and  X where  X and 

X include  all  nodes.  When  two  nodes  are  specified  (such  as  an  input 

N-2 

and  an  output  node)  the  number  of  proper  cuts  is  reduced  to  2 

Jensen  (2)  presents  an  algorithm  that  uses  a tree  search  approach  for 

determining  proper  cuts.  The  generation  of  proper  cuts  is  a lengthy 

one  and  Jensen  proposes  tlje  use  of  restricted  proper  cuts  to  decompose 

the  graph.  A restricted  proper  cut  is  defined  as  a set  of  nodes  X such 

that  if  there  is  an  arc  from  node  i to  node  j , and  if  node  j belongs  to 

X,  then  node  i must  also  belong  to  X. 

The  network  of  Figure  1 has  187  proper  cuts  and  69  restricted  proper 

13 

cuts.  A network  with  thirteen  nodes  could  have  as  many  as  2 or  8192 

proper  cuts.  The  network  can  be  divided  into  1533  modules  by  using  proper 

cuts  and  only  696  modules  by  using  restricted  proper  cuts.  There  could 
13 

be  as  many  as  E ^ S(13,n)  cuts  where  S(13,n)  is  the  Sterling  number 

of  the  second  kind  (S(m,n)  * E^^-l)1"^)  (n-i)m/n!). 

Suppose  the  network  in  Figure  4a  must  be  partitioned.  Nodes  1,  2, 
and  3 cannot  be  in  a single  module;  however,  any  pair  of  nodes  can  be  in  the 
same  module.  Let  node  1 be  the  source  node  and  node  3 be  the  sink  node. 

The  penalty  for  cutting  arc  (1,3)  is  very  high,  say  M,  whereas  the  penalty 
for  cutting  arcs  (1,2)  and  (2,3)  is  very  small,  say  m,  where  2m  « M. 
Decomposition  by  proper  cuts  would  force  arc  (1,3)  to  be  cut  and,  there- 
fore, experience  the  high  penalty  (M+m).  The  network  would  be  optimally 
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decomposed  by  making  two  sets  of  nodes,  set  1 contains  node  1 and  node  3 
and  set  2 contains  node  2 at  a small  penalty  (2m)  , which  is  very  much 
less  than  the  penalty  for  cutting  arc  (1,3).  This  partitioning  by  proper 

cuts  is  shown  in  Figures  4a  and  4b.  The  optimal  partitioning  is  shown 
in  Figure  4c. 

V.  LIFE  CYCLE  COST  MODEL 

Life  Cycle  Cost  (LCC),  as  used  in  this  work,  is  the  cost  of  acquisition 
plus  the  suport  cost  for  an  equipment  over  its  intended  useful  life.  The 
LCC  model  used  is  after  Caponecchi  (4)  and  is: 

LCC  * Ct(M1)  + Ct(M2)  + ...  + Ct(Mi)  +...+  Ct(Mn), 
where  LCC  *=  Life  Cycle  Cost 

= Module  i of  design  m 

CtfM^  = The  cost  of  acquiring  and  maintaining 
module  i generated  by  design  m. 
n = Total  number  of  modules  in  design  m. 
or  LCC  - Z1"1  CtO^).  (1) 

Since  the  Life  Cycle  Cost  is  composed  of  two  major  elements;  namely,  the 
acquisition  cost  and  the  life  time  support  cost;  the  acquisition  cost  (not 

including  spares)  can  be  expressed  in  terms  of  the  candidate  modules.  That 
expression  is: 

C*  - »*  E1?1  C(Mt>  (2) 

Where  Ca  «=  Acquisition  Cost  for  Ne  equipments 
Ne  = Number  of  equipments  to  be  procured 
C(Mt)  - Cost  of  acquiring  module  i of  design  m. 
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Under  the  assumptions  of  a Discard  at  Failure  Maintenance  (DAFM) 

policy  and  that  at  least  one  of  each  module  will  be  spared,  the  life 

time  support  cost  of  Ne  equipments  if  each  is  at  a separate  site  is: 

Ce  uCi  + nCcL  + Ne  E " N C(M.)  (3) 

i=l  1 

Where  Cs  = Total  Organization  Support  Cost 

Ci  = Cost  of  introducing  a line  item  into  the 

5 . 

inventory  system 

Cc  = Cost  of  maintaining  a line  item  in  inventory 
for  one  year 

L = Planned  operational  life  of  the  equipment, 
in  years 

= Number  of  spares  of  module  i to  be  procured 
to  support  each  equipment 

The  total  Life  Cycle  Cost  is  the  sum  of  Equations  (2)  and  (3),  or: 

LCC  = n(Ci  + CcL)  + Ne  E^  (1+1^)  CO^)  (4) 

The  Life  Cycle  Cost  for  a single  equipment  is 

LCCe  = N1(Ci+CcL)  + E±”x  (1+N^  C(Mi>  (5) 

Where  N^  = n/Ne. 

Equation  (5)  assumes  that  the  two  major  costs  that  are  affected  by 
the  modularization  of  a design  are  the  acquisition  cost  and  the  organ- 
ization support  cost  associated  with  the  physical  inventory  of  parts. 

It  ignores  the  cost  of  maintenance  personnel  and  consumables  on  the 
assumption  that  they  will  be  relatively  unaffected  by  the  modularization 
decisions.  It  further  assumes  there  will  be  "logistic  self-sufficiency". 
The  implication  being  that  the  shelf  life  of  all  items  in  the  spares 
inventory  exceeds  the  planned  life  of  the  equipment  and  that  there  is  no 
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deterioration  during  non-use.  A further  assumption  is  that  there  will 
be  only  one  level  of  modularization. 

NOTE:  If  Equation  (4)  is  minimized  on  the  basis  of  an  integer  number 
of  spares  for  each  equipment,  is  integer.  If  the  number  of 
spares  for  all  equipments  is  to  be  integer,  then  NeN^  must  be 
Integer. 


VI.  ADEQUACY  OF  SPARES  CONSTRAINT 

All  modules  are  assumed  to  have  a constant  failure  rate  while  in 
service  and  a zero  failure  rate  while  in  spares.  With  these  assumptions, 
the  probability  that  N^  spares  will  be  adequate  to  meet  the  availability 
over  the  planned  life,  L,  of  the  equipment  is: 

-f, 


where 


P(w±  < N±) 


if  Wi 
e f . 

E nn  r > A 

w“0  w^.  — 


P(w±  N^)  * Probability  that  the  number  of  failures 
of  module  i is  no  larger  than  the  number 
of  spares  of  module  i. 
w.  * Number  of  failures  of  module  1 in  L. 


(6) 


f ^ « Expected  number  of  failures  of  module  i in  L. 

A - The  minimum  allowable  availability  of  an 
equipment . 


An  expression  for  f^  is: 


f±  - rt  L Ne, 


(7) 


where  r^  “ failure  rate  of  module  1. 

The  failure  rate  of  module  1,  r^,  is  sensitive  to  the  modularization  since 
the  failure  rate  of  a module  is  a function  of  the  failure  rates  of  the 
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elements  in  the  module  and  the  number  of  external  and  internal  connections. 
This  can  be  expressed  as: 


jeM.  j 


r . - AP  (M. ) r , 


i P’ 


where  r^  = The  failure  rate  of  element  j of  module  i 

AP^)  = The  number  of  external  connections  eliminated  by 

putting  elements  j in  module  i. 

r = Failure  rate  correction  factor  for  interconnection 
P 

reduction. 

An  important  point  to  be  noted  is  that  the  value  of  r^  in  Equation  (8) 
is  linearly  dependent  on  the  number  of  external  connections  eliminated 
by  the  modularization. 


VII.  PHYSICAL  CONSTRAINTS 

In  the  modularization  of  equipments,  there  will  be  some  physical 
constraints  which  cannot  be  violated.  In  electronics  modularization, 
the  constraints  might  be:  weight,  heat  generated,  external  connections, 
number  of  elements,  chip  area,  power  consumption,  physical  and  electronic 
compatibility,  etc.  The  exact  nature  of  the  constraints  is  not  as 
important  as  their  form.  For  example,  linear  constraints  are  usually 
much  easier  to  handle  than  non-linear  constraints.  However,  non-linear 
constraints  that  are  convex  will  usually  give  less  trouble  than  non-linear 
constraints  that  are  neither  concave  nor  convex.  Caponecchi  (A)  has 
assumed  some  linear  constraints  which  do  not  cause  any  difficulty.  Such 
physical  constraints  as  are  mentioned  above  can  be  expressed  as: 


VMj  V-Aik"“ 


1!) 


where 


- — — 


* The  value  of  characteristic  k for  element  j 
assigned  to  module  i. 

A^max  ■ Maximum  value  of  characteristic  k that  may  be 
assigned  to  module  1. 

VIII.  MAINTENANCE  CONSTRAINT 

The  expected  maintenance  time  for  a module  will  be  a function  of 
the  size  of  the  module  (the  number  of  elements  (SFE’s)  in  the  module) 
and  the  number  of  external  connections  to  the  module.  Caponecchi  (4) 
proposed  that  the  maintenance  time  is  composed  of  the  sum  of  a constant 
time  element,  a time  element  that  is  a linear  function  of  the  number  of 
modules  in  the  equipment  and  a third  time  element  that  is  exponentially 
related  to  the  number  of  external  connections  to  the  module.  A general 
expression  for  the  expected  maintenance  time  is: 

E(TM)  - + nT2  + E1"1  T3  expC^P^))  (10) 

where  E(TM)  “ Expected  maintenance  time. 

T^  ■ A constant  time  per  maintenance  action. 

T^  ■ A constant  time  per  module. 

T^  ■ A constant  modifying  the  exponential  relationship 
of  the  number  of  external  connections,  and 

T^  ■ A constant  modifying  the  number  of  external  connections. 

P(M^)  • Number  of  external  connections  to  module  i. 

The  specific  form  of  Equation  (10)  used  by  Caponecchi  (4)  is: 

E(TM)  - 2.5  + 0.05n  + 0.087  Z ^ exp(0.047  P(M1)).  (11) 

This  expression  yields  a constraint  on  the  expected  time  for  maintenance 
which  must  be  less  than  the  maximum  mean  time  to  repair;  or 

E(TM)  < MTTRmax.  (12) 

where  MTTRmax  - maximum  allowable  mean  time  to  repair. 
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IX.  MATHEMATICAL  FORMULATION  OF  THE  MODULARIZATION  PROBLEM 


From  the  previous  discussion,  it  is  now  possible  to  formulate  the 
modularization  problem.  That  formulation  over  all  equipments  is: 

Min.  LCC  - n(Ci+CcL)  + Ne  (1+N±)  CCM^ 


Subject  to:  -f^ 


2) 

3) 


T1  + nT2  + T3  Ei-1  exP <T4P (**!>)  i MTTRmax 

E.  u A..  < A.  . max  for  all  i and  k 
jeM^  jk  — i-k 

> 0 and  Integer. 


(13) 


The  formulation  for  each  equipment  is: 

Min  LCC  - N^Ci+CcL)  + Z (1+N±)  C(M  ).  (14) 

This  expression  for  the  LCC  is  subject  to  the  same  constraints  as  listed 
above  for  Equation  (13). 


X.  DEPENDENCE  OF  LCC  ON  EXTERNAL  CONNECTIONS 

The  mathematical  formulation  of  the  modularization  problem  shows 
that  the  solution  is  dependent  on  the  number  of  external  connections  in 
three  ways.  First,  the  minimum  LCC  (objective  function)  is  dependent 
on  the  number  of  external  connections  through  C(M^),  the  cost  of  a module. 
Since  internal  connections  are  less  expensive  than  external  connections, 
the  cost  of  modularizing  N elements  into  modules  will  be  reduced  as  some 
function  of  the  reduction  of  external  connections.  The  cost  of  packaging 
is  also  reduced.  Second,  the  failure  rate  of  a module  is  a function  of 
the  number  of  external  connections  as  shown  by  Equation  (8) . The  avail- 
ability constraint  is  a function  of  the  failure  rate  of  the  modules.  And, 
third,  the  maintenance  constraint  is  a function  of  the  number  of  external 


connections  through  its  exponential  term.  Others  have  found  the  LCC 
to  be  highly  depeAlent  on  the  reliability.  Hardy  states  "...Since  hard- 
ware failure  frequency  is  one  of  the  key  parameters  that  influences  life- 
cycle  cost..."* 

A further  consideration  is  that  there  is  no  guarantee  that  there 
exists  a modular  arrangement  that  satisfies  all  constraints.  Since  the 
physical  constraints  restrict  the  size  of  modules  and  availability  will 
increase  with  the  size  of  the  modules,  the  availability  constraint  may 
be  set  at  a value  that  is  unattainable  in  any  modularization  process.  .'he 
same  may  be  true  for  the  mean-time-to-repair  constraint. 

XI.  A MODIFIED  PROBLEM 

Since,  as  was  shown  in  the  previous  section,  the  Life  Cycle  Cost  is 
a function  of  the  number  of  external  connections  between  modules,  the 
partitioning  or  modularization  method  will  be  based  on  the  reduction  in 
external  connections  achieved  by  the  modularization.  Assume  the  expected 
cost  of  an  external  connection  is  Ce  and  the  expected  cost  of  an  Internal 
connection  is  Ci  where  Ce  is  greater  than  Ci.  The  modified  objective 
function  is  to  minimize  the  cost  for  connections;  that  is: 

Min  C * (E1Ce  + I Ci)  (15) 


*Hardy,  C.A.  "Avionic  Reliability  and  Life-Cycle-Cost  Relationship", 

AGARD  Lecture  Series  No.  81  on  Avionics  Design  for  Reliability. 

(Available  from  NASA,  Langley  Field,  Va.  23365,  Attn:  Report  Distribution 
and  Storage  Unit.) 


Where  C * Cost  of  connections  to  and  within  a module, 

= Number  of  external  connections  from  module  i, 

Ce  = Cost  of  an  external  connection  to  a module, 

1^  = Number  of  internal  connections  in  module  i,  and 
Ci  ■ Cost  of  an  internal  connection  in  a module. 

The  constraints  on  the  physical  design  parameters  will  remain  the  same. 

The  availability  of  the  equipment  is  increased  as  the  reliability  is 
increased  and  the  mean-time-to-repair  is  decreased  as  the  number  of 
modules  and  external  connections  is  decreased.  Thus,  if  the  size  of  a 
module  is  at  its  maximum  value,  its  availability  will  tend  toward  a 
maximum  and  its  MTTR  will  tend  toward  a minimum. 

The  modularization  process  proceeds  to  combine  SFE's  in  a sequential 
manner  until  no  other  SFE  can  be  added  without  violating  a constraint. 

The  module  so  constructed  is  set  aside  and  a new  module  constructed  from 
the  remaining  SFE's.  The  process  continues  until  all  SFE's  have  been 
assigned  to  a module.  Since  the  assignment  of  elements  to  modules  under 
this  procedure  will  be  dependent  on  the  starting  point,  it  was  decided 
to  do  N iterations  with  one  iteration  starting  from  each  node.  The  best 
modular  arrangement  is  then  chosen  from  the  N arrangements.  (It  is 
expected  that  the  N arrangements  are  not  unique  so  that  fewer  than  N 
modular  arrangements  need  to  be  evaluated.)  The  best  module  being  the 
one  with  the  smallest  difference  between  external  connections  and  internal 
connections.  (This  number  can  be  negative.) 

Once  the  best  module  is  chosen,  it  must  be  evaluated  for  its  mean- 
time-to-repair.  If  it  does  not  meet  the  mean-time-to-repair  constraint, 
the  second  best  module  is  evaluated,  then  the  third  best,  etc.  until  a 
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. Put  Mj  aside,  put  j * j + 1. 

9.  Remove  the  node  most  recently  added  to  from  KN  '.  If  MN 
Is  null,  go  to  Step  10.  Otherwise  go  to  Step  5. 

1C.  After  all  nodes  have  been  assigned  to  modules,  store  the  modules 
as  set  Sk.  Put  k * k+1.  Put  1 = 1+1.  If  1+1  > N,  go  to  Step  11. 
Otherwise  go  to  Step  1. 

11.  Find  that  set  of  modules  with  the  smallest  value  of  value  of  ST 
when  summed  over  all  modules  In  the  set. 

12.  Calculate  the  MTTR.  If  the  MTTR  <_  MTTRmax,  STOP.  This  is  the 
solution.  Otherwise  remove  current  set  of  modules  from  further 
consideration  and  go  to  Step  11.  If  no  set  of  modules  satisfies 
the  MTTR  constraint,  design  Is  infeasible. 

XIII.  CALCULATION  OF  THE  NUMBER  OF  EXTERNAL  AND  INTERNAL 
CONNECTIONS  FOR  THE  MODULES 

A network  can  be  represented  by  a matrix  and  this  is  the  representation 
used  to  calculate  the  values  of  EXT,  INT  and  ST  in  previous  solution  pro- 
cedure. The  matrix  representation  of  the  network  in  Figure  1 (without  input 
and  output)  is  shown  in  the  matrix  in  Figure  5.  The  entries  in  row  i are 
the  connections  out  of  i and  entries  in  column  j are  the  connections  into  j. 
The  entries  on  the  main  diagonal  are  the  total  number  of  connections  to 
node  1 (both  into  and  out  of  node  i). 

The  determination  of  the  best  node  to  associate  with  node  1 is  shown 


in  Table  1 


Figure  5: 
Table 

Node 

Pair 

L,2 

L.3 

L,4 

L.5 

L,6 

M 

L»8 
L,9 
L,  10 

Ml 
M2 
M3 


Matrix  Representation  of  Network  Shown  in  Figure  1 
1.  Determination  of  Best  Node  to  Associate  With  Node  1 
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22 

32 
14 
22 
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Table  1 can  be  found  as  follows: 


EXT  = Sum  of  entries  in  row  1 (except  for  main  diagonal  element) 

+ sum  of  entries  in  row  2 (except  for  main  diagonal  elements) 
- the  entry  at  1,  2. 

* (5+6+1)  + (3+2+3)  - 5 = 15 
INT  = Entry  at  intersection  1,  2 
■=  5 

ST  - EXT  - INT  - 15  - 5 = 10 

The  minimum  difference  between  the  number  of  external  connections  and 
the  number  of  internal  connections  occurs  when  nodes  1 and  3 are  paired. 

The  pairing  of  nodes  1 and  3 must  be  checked  against  the  constraints. 

The  physical  characteristics  of  the  nodes,  as  used  by  Caponecchi,  is 
given  in  Table  2. 

The  combination  of  nodes  1 and  3 has:  59  Part  Count,  385  Heat  Generation, 

.049  Required  Chip  Area  and  9 External  Connections.  The  maximum  allowable 
values  are  150,1100,0.17  and  25  respectively.  So  Nodes  1 and  3 should  be 
combined  into  a module. 

The  matrix  can  be  reconstructed  by  combining  nodp»  1 and  3 into  node 
1'  . Two  entries  are  made  in  the  1,'  1 'location  the  values  are  the  total 
number  of  external  connections/the  number  of  internal  connections.  The 
matrix  with  this  entry  is  shown  in  Figure  6.  Using  Figure  6,  the  best 
node  to  be  combined  with  1,  3 is  2.  The  value  of  ST  for  nodes  1,  3,  2 
in  a module  is  1 (the  external  connections  are  (3+2+3+3+1)  and  the 
internal  connections  are  (5+6)).  The  computations  are  shown  in  Table  3. 


Table  2.  Physical  Characteristics  of  the  Nodes* 


■ 


■ 


i 


Node 

Part 

Count 

Heat 

Gneration 

Required 
Chip  Area 

Number  of 

External 

Connections 

Failure  Rate/Hr. 
x 10“  ^ 

1 

27 

180 

.024 

12 

.211 

, 2 

43 

210 

.026 

13 

.450 

3 

32 

205 

.025 

9 

.225 

4 

40 

290 

.038 

5 

.300 

5 

25 

220 

.020 

13 

.190 

6 

15 

140 

.022 

10 

.178 

7 

27 

185 

.024 

20 

.420 

8 

18 

150 

.020 

2 

.080 

9 

23 

200 

.022 

10 

.160 

10 

50 

260 

.050 

22 

.350 

11 

26 

150 

.023 

9 

.182 

12 

17 

145 

.020 

8 

.140 

13 

12 

80 

.015 

10 

.170 

♦maximum 
allowable 
per  module  150 

1100 

0.17 

25 

19 


•;a 

• v 


-.•3 

i • 


First  Reduced  Matrix 


Table  3.  Determination  of  Best  Node  to  Associate  With 
The  Node  1,3  Combination 


1 Minimum  ST 


The  physical  characteristics  are: 

Part  Count  102  < 150 
Heat  Generation  595  < 1100 
Required  Chip  Area  .075  < .17 
Number  of  External  Leads  12  < 25. 

So  a module  can  be  made  from  nodes  1,  2 and  3. 

Proceeding  in  the  above  manner,  the  best  modules  starting  from 
node  1 are: 

- (1,2,3)  (EXT=15,  INT=12) 

M2  = (5,6,9,10,12,13)  (EXT=  14,  INT=29)  and 
Kj  - (7,8,11),  (EXT=19 , INT=6) 

The  total  number  of  external  connections  is  48  and  the  number  of  internal 
connections  is  47.  The  value  of  ST  is  1. 

The  solution  process  now  performs  an  additional  N-l  iterations 
yielding  N designs.  These  designs  will  probably  not  be  unique,  in  fact, 
for  the  example  problem,  there  are  only  5 unique  modularizations.  These 
five  designs  are  shown  in  Table  4.  It  should  be  noted  that  there  are  two 
designs  with  the  same  minimum  number  of  external  connections  and  that  those 
modular  arrangements  are  obtained  by  starting  from  any  of  the  first  8 nodes. 
Under  the  assumed  model,  the  mean- 1 ime- to-r epair  will  be  minimum  when  the 
number  of  external  connections  is  minimum.  It  now  becomes  necessary  to  see 
if  the  MTTR  of  the  designs  created  by  the  modularization  process  meet 
the  MTTRmax  requirement. 


Table  4.  Five  Unique  Designs 


Starting  Nodes 

Elements  in 
Modules 

Total  // 

External  Connections 

1,2,3 

1,2, 3, 4 

5,6,9,10,12,13 

7,8,11 

48 

4, 5, 6, 7, 8 

3,4,7 ,8,11 
5,6,9,10,12,13 

1,2 

48 

9 

5,8,9,11,12,13 

4,6,7,10 

1,2,3 

52 

10 

2,5,6,10 

4,8,9,11,12,13 

1,3,7 

58 

11  12  13 

5,8,9,11,12,13 

XIV.  CHECKING  MTTR  CONSTRAINTS 

The  model  assumed  that  the  expected  time  for  maintenance  took  the 

form: 

E(TM)  = 2.5  + 0. 05n  + 0.087  Z ^ expC.OATPCM^^)) . (11) 

The  MTTR  will  be  checked  for  conformance  to  the  MTTR  constraint.  The 


MTTR  has  been  set 

max 

at  3.36. 

The 

calculations 

are  shown  in  Table  5. 

Table  5. 

MTTR 

Calculations 

Module 

EXT 

INT 

ST 

exp(.047P(M1)) 

(1,2, 3, 4) 

15 

12 

3 

2.0239 

(5,6,9,10,12,13) 

14 

29 

-15 

1.9309 

(7,8,11) 

19 

6 

13 

2.4572 

E(TM)  = 2.5  + 0.15  + 0.558  = 3.208  < 3.36. 

Thus,  from  Table  5,  it  can  be  seen  that  the  modularization  does  meet 
the  MTTR  constraint. 


XV.  SPARES  REQUIREMENT  TO  SATISFY  AVAILABILITY 

Since  every  element  must  function  for  the  system  to  function,  the 
modules  can  be  treated  as  a series  system  from  the  reliability  stand- 
point. The  availability  constraint  is  dependent  on  the  module  reliability 
in  accordance  with  Equations  (6),  (7),  and  (8).  That  is: 


— f v 

i * i 
e f . 

pfw  < N ) * £ ° „ ? > A 

rtwi  _ w *0  w±!  - 


f^  - r^  L Ne,  and 


ri  - ^ rj  -*P(VV 


(6) 

(7) 

(8) 


Where  r - 1 X 10 
P 


-6 
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The  cost  of  a module,  cOl^),  is: 

C(M.)  = C,  N + C_  EXT  + C, 

1 l p i.  i 

Where  “ Components  parts  cost  factor 

N “ Number  of  parts  in  the  module 
P 

0,^  = Interconnection  cost  factor 
EXT  ■ Number  of  interconnections 
= Packaging  cost 

The  failure  and  cost  data  are  summarized  in  Table  6.  These  data  are 
for  the  design  that  starts  at  nodes  1,  2,  or  3.  The  total  cost  for 


design  is  $1,106. 

This  design  has 

the  least  cost  of 

any  of  the 

Table  6. 

Failure  and  Cost 

Data  For  a Modular  Design 

Module 

Failure  Rate 
106  Hrs. 

Expected  Number 
of  Failures 

Cost  of 
Module 

(1.2, 3,4) 

26.9 

71.0 

$ 414 

(5,6,9,10,12,13) 

25.9 

68.0 

412 

(7,8,11) 

7.8 

20.5 

280 

designs  shown  in  Table  4. 

The  availability  problem  can  be  expressed  as: 

Min  C - 414NX  + 412N2  + 280N3 

subject  to: 

N1  * N2  x 

(E^J  (exp(-71))(71x)/x!)*(E^0  (exp(-68)) (68x)/x!) . 

N3  x 

(E^0(exp(37))(37x)/x:  > 85  (16) 

N^N^Nj  0 and  Integer 


J 


Since  the  above  problem  is  an  integer  non-linear  programming  problem. 


an  approximate  solution  method  will  be  used.  The  solution  method  is: 

1)  Initialize  by  finding  a lower  bound,  , for  each  module. 
Choose  Nj  such  that 


N 

Zx^0  (exp(~fi))(fix)/x!  > A = .85 

2)  Find  the  value  of  A for  the  system. 


Where  A *=  ni”1(Ix^0(exp(-fi))  (f ix) /x! 

If  A >_  .85,  STOP. 

3)  Let  N±  = N±  + 1. 

N +1 

4)  Calculate  AA.  = (exp(-fi))f.x/x!)  - 

<£„  W-f,)!,’/-!. 

5)  Calculate  AA^/C^) 

6)  Choose  max  AA  /C(M  ).  Set  all  other 

i 1 1 


- 1,  Co  to  Step  2. 

The  calculations  for  the  example  are: 

Nt  - 80,  N?  - 77,  N3  - 43. 

And  A = (.8692)(. 8741) (.8568)  - .6509.  Now  set  N - 81,  N£  - 78  and 
N3  - 44.  The  values  of  AAj  are:  AA^  - .^226,  AA^  - .0224,  AA3  * .0321 
and  the  values  of  AA^CtOy  are:  AAJ/Ct(M1)  - .0000545,  AA2/Ct(M2)  - 
.0000543  and  AA3/Ct(M3)  - .0001146.  So,  with  N3  - 44  and  Nj  = 80  and 
with  N2  * 77,  calculate  A (A  ■ .6753).  Since  A is  less  than  0.85  it 
is  necessary  to  perform  other  iterations  until  A > .85.  The  final 


i 
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values  for  are  * 84,  ^ = 80,  and  = 49.  This  process  requires 
thirteen  iterations  and  the  final  value  of  A is  A * 0.8572. 

The  LCC  for  the  system  depends  on  the  spares  requirements  that  are 
necessary  to  meet  the  availability  requirement.  However,  the  minimization 
of  the  LCC  is  independent  of  the  determination  of  the  spares. 

XVI.  COMPUTER  PROGRAM 

A computer  program  has  been  written  to  accomplish  the  above  described 
modularization.  An  explanation  of  the  program  and  its  flow  chart  are  pre- 
sented as  Appendix  1.  The  program  listing  is  presented  as  Appendix  2. 

Some  output  for  the  example  problem  is  included  as  Appendix  3. 

XVII.  SUMMARY 

This  research  has  developed  a method  for  modularizing  an  equipment 
subject  to  physical  constraints,  a mean-time-to-repair  constraint  and 
an  availability  constraint.  The  model  for  life  cycle  cost  that  was  used 
is  highly  dependent  on  the  physical  configuration;  namely,  the  number  of 
external  connections  to  a module. 

Although  the  modular  arrangements  developed  have  not  been  proven  to 
be  optimal,  no  counter-examples  have  been  seen.  Further  work  would  be 
desirable  to  1)  prove  the  method  proposed  is  optimal,  and  2)  if  1)  cannot 
be  achieved,  then  develop  an  optimal  mathematical  solution  to  the  problem. 

The  solution  method  and  assumed  model  do  not  require  a determination 
of  spares  requirements  while  performing  the  modularization  process.  It  is 
recognized,  though,  that  the  cost  of  a module  could  be  dependent  on  the 
total  number  of  items  purchased  and  in  such  a cases  spares  would  enter 
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into  the  minimization  process.  The  necessity  for  buffers  between 
stages  and  a strong  desire  to  group  certain  components  into  a given 
module  were  not  addressed. 
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APPENDIX  I:  Description  of  the  Computer  Program 

The  network  partitioning  program  resides  in  the  computer  in 
the  Multics  environment.  The  program  including  its  subroutines  and 
the  network  data  are  stored  in  segments.  .his  enables  the  user  to 
access  any  part  of  the  program  for  editing  or  compilation  by  calling 
the  subroutine  by  name. 

The  names  for  the  segments  that  contain  the  program  parts  are: 

1)  netpart. fortran 

2)  minx. fortran 

3)  sortx. fortran 

4)  sortxl . fortran 

5)  mshift. fortran 

6)  rshift. fortran 

7)  network. data 

A brief  description  of  the  subroutines  follows: 

1)  netpart--contains  the  coding  for  the  algorithm  described  in 
Section  XII  (pages  15  and  16).  (The  program  currently  re- 
siding in  the  RADC  computer  can  partition  a network  that  has 

as  many  as  eighty  nodes  and  as  many  as  ten  constraints.) 

2) minx--contains  a subroutine  called  by  netpart  that  selects 

the  group  of  elements  with  the  minimum  value  for  the  external 
minus  the  internal  connections  at  each  iteration. 

3) sortx  and  sortxl--contain  subroutines  to  properly  arrange  the 
identification  of  the  elements  in  each  module. 

4) mshif t--contains  a subroutine  to  rearrange  the  data  matrix  so 
that  every  node  is  considered  as  a starting  node  for  a set  of 
modules . 

5) rshif t--contains  a subroutine  to  reassign  physical  property 
values  after  the  matrix  has  been  rearranged  by  mshift. 

6) network.data--contains  all  of  the  data  relative  to  a given 
network. 
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It  is  assumed  that  the  user  can  create  segments  where  the  data  on 
the  network  to  be  partitioned  can  be  stored.  Data  for  different  networks 
can  be  stored  in  different  segments  with  different  names.  The  data 
for  the  example  network  is  stored  in  a segment  called  ' network. data ' . 

When  a user  wants  to  store  data  it  is  necessary  to  enter: 

'edm.  "network  data  name"  ' 

If  the  segment  name  already  exists,  the  used  will  be  put  into  the  edit 
mode.  If  a new  segment  is  being  defined,  the  user  will  be  put  into  the 
input  mode  and  can  enter  the  data  on  the  new  network.  Fig.  1-1  shows 
the  procedure  for  creating  a 'data  segment.  An  explanation  follows: 

Line  1--  standard  login  procedure. 

Line  10--user  calls  an  already  existing  segment  so  Multics 
puts  the  system  into  the  edit  mode. 

Line  11 — * Q 1 is  the  quit  command.  No  change  is  made  in  the  data 

and  Multics  acknowledges  by  giving  the  system  status  data. 

Line  13--user  calls  the  edm  function  with  a new  name  1 edm . ex. da ta ' 
Multics  searches  far  such  a segment  andsincfc  it  doesn't 
find  such  a segment,  responds  with  the  message  on  line  14. 

Line  15--indicates  to  the  user  to  input  his  data. 

Line  16--type ' . ' to  change  from  input  to  edit  mode  after  data  .has 
been  entered. 

Line  17--Multics  response  to  the  request  to  enter  the  edit  mode. 

Line  18--type  'Q'  to  quit  without  storing  data,  otherwise  type 
'W'  with  a carriage  return  and  ' Q ' to  store  new  data. 

The  printout  of  the  data  for  the  example  network  is  shown  in  Fig  1-2. 

The  input  data  should  be  put  in  in  this  same  format.  An  explanation 
of  Fig.  1-2  follows. 

Line  1--  number  of  nodes  in  the  network,  N. 

Line  2--  number  of  physical  constraints  on  the  partitioning,  k. 

Lines  3 ,4&5- -vec tors  Cl  where  1=1, 2, ...,k.  Each  vector  has  N 
elements  and  there  are  k vectors* 

Lines  6 through  18--the  number  of  connections  between  nodes  or 
elements . 

Line  19--a  vector  that  contains  the  limits  on  the  vectors,  Cl. 
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After  the  user  has  entered  the  data  for  the  network  and  attached 
it  to  the  main  program,  the  program  is  ready  to  be  compiled  and  run. 
The  output  gives  the  best  modular  arrangement  starting  from  each  of 
the  N nodes.  The  output  also  includes  the  value  of  ST  (external  minus 
internal  connections)  for  each  module.  The  user  can  then  compare  the 
N arrangements  to  determine  the  'best'  configuration.  The  output  for 
a sample  problem  is  shown  in  Appendix  III. 


Fig.  1-1.  Creating  a Data  Segment 


Printout  of  a Data  Segment 


APPENDIX  II:  Flow  Chart  for  the  Computer  Program 

Appendix  II  contains  the  flow  chart  for  the  computer  program 
previously  described.  The  flow  chart  is  presented  in  Fig.  II— 1 - 


Fig.  1 1 - 1 : Flow  Chart  for  Computer  Program 
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APPENDIX  III:  An  Example  of  the  Output  of  the  Computer  Program 

The  output  from  the  computer  program  for  the  example  problem 
is  shown  in  Fig.  I I I— 1 . An  explanation  of  Fig.  III-L  follows. 

Lines  1-9- -standard  login  procedure 

Line  10 — asks  the  Fortran  compiler  to  compile  the  main  program 
segment  containing  netpart. 

Line  13--the  user  asks  the  compilation  of  all  of  the  subroutines 
by  typing  the  name  of  each  segment.  Lines  16,19,  22  and 
25  compile  all  of  the  segments. 

Line  28--the  input/output  command  which  attaches  the  data  to  the 
main  program. 

Line  30--requests  output. 

Line  32--starts  the  computation. 
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METRIC  SYSTEM 


BASE  UNITS: 

Quantity 


SI  PREFIXES: 


Multiplication  Factors 


1 000  000  000  000  - 
1 000  000  000  = 
1 000  000  « 
1 000  - 
100- 
10- 
0.1  = 
0 01  = 
0.001  - 
0.000  001  - 
0 000  000  001  - 
0.000  000  000  001  - 
0.000  000  000  poo  001  - 
0.000  000  000  000  000  001  - 

i avoided  where  possible 


SI  Symbol 


Formula 


length 

metre 

m 

mass 

kilogram 

kg 

time 

second 

s 

electric  current 

ampere 

A 

thermodynamic  temperature 

kelvin 

K 

amount  of  substance 

mole 

mol 

luminous  intensity 

candela 

cd 

SUPPLEMENTARY  UNITS: 

plane  angle 

radian 

rad 

solid  angle 

steradian 

sr 

DERIVED  UNITS: 

Acceleration 

metre  per  second  squared 

m/s 

activity  (of  a radioactive  source) 

disintegration  per  second 

(disintegration)fs 

angular  acceleration 

radian  per  second  squared 

rad/s 

angular  velocity 

radian  per  second 

rad/s 

area 

square  metre 

m 

density 

kilogram  per  cubic  metre 

kg/m 

electric  capacitance 

farad 

F 

A-s/V 

electrical  conductance 

siemens 

S 

A/V 

electric  field  strength 

volt  per  metre 

V/m 

electric  inductance 

henry 

H 

V-s/A 

electric  potential  difference 

volt 

V 

W/A 

electric  resistance 

ohm 

V/A 

electromotive  force 

volt 

V 

W/A 

energy 

joule 

1 

N-m 

entropy 

joule  per  kelvin 

)/K 

force 

newton 

N 

kg-m/s 

frequency 

hertz 

Hz 

(cycle)/s 

illuminance 

lux 

lx 

Im/m 

luminance 

candela  per  square  metre 

cd/m 

luminous  flux 

lumen 

Im 

cd-sr 

magnetic  field  strength 

ampere  per  metre 

Wb 

A/m 

magnetic  flux 

weber 

V-s 

magnetic  flux  density 

tesla 

T 

Wb/m 

magnetomotive  force 

ampere 

A 

)/* 

power 

watt 

W 

pressure 

pascal 

Pa 

N<m 

quantity  of  electricity 

coulomb 

C 

A-s 

quantity  of  heat 

joule 

1 

N-m 

radiant  intensity 

watt  per  steradian 

W/sr 

specific  heat 

joule  per  kilogram-kelvin 

I'kg-K 

stress 

pascal 

Pa 

N/m 

thermal  conductivity 

watt  per  metre-kelvin 

W/m-K 

velocity 

metre  per  second 

m/s 

viscosity,  dynamic 

pascal-second 

Pa-s 

viscosity,  kinematic 

square  metre  per  second 

m/s 

voltage 

volt 

V 

W/A 

volume 

cubic  metre 

m 

wavenumber 

reciprocal  metre 

(waveym 

work 

joule 

1 

N-m 
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MISSION 

of 

Rome  Air  Development  Center 


RADC  plans  and  conducts  research,  exploratory  and  advanced 
development  programs  in  commeLnd,  control,  and  coimunications 
(C3)  activities,  and  in  the  C3  areas  of  information  sciences 
and  intelligence.  The  principal  technical  mission  areas 
are  communications,  electromagnetic  guidance  and  control, 
surveillance  of  ground  and  aerospace  objects , intelligence 
data  collection  and  handling,  information  system  technology, 
ionospheric  propagation,  solid  state  sciences,  microwave 
physics  and  electronic  reliability , maintainability  and 
compatibility . 
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